Visual Basic program "VB 2019 Sagittarius"

VB2019 Simultion of our Galagy, including sagitarius and 10 mini Black holes

Introduction and Purpose

The purpose of the program "VB Sagittarius" is to simulate the black hole Sagittarius A* and 10 mini BH's:
See also:

This program is part of what is called the "Sagitarrius Project". To read more select this link Sagittarius Project
Temporary important is this: Carpe Diem and Sagittarius.
Please study that document in relation with Display 2, shown below, and observe that the w parameter changes when Display 2 is selected or not selected. The same is true for inclination. The omega parameter is constant, which it should not.

Program Description

The program consists of 4 Forms: Control, Display, Menu and Figure

The Control Form is used to select a planet simulation and to perform a simulation.
The Display Form shows the results of a simultion and the result of a simulation saved in a Data Bae
The Menu Form is to select additional parameters used by a planet simulation
The Figure Form is to select a simulation saved in a Data Base

Control Form - Operation

Display 1

The control form is used to select and perform a simulation.

First select a particular simulation to be performed.
In the example this is simulation 7, which demonstrates "Sagittarius and 10 Mini BH's"
Next Select "Parameter Selection Menu" in order to change cetain parameters relevant for this simulation. This is an option.
Next Select Start
Next the simulation will show the Display Form as shown in Display 2. An additinal display shows the evolution of the simulation.See Display 3.

The control form is also use to create a certain number of predefined figures.
This option is currently not available.

First select a particular figure to be created.
Select Menu Selection Form. Select SubTest and return back to control form
Next Select Start
Next the simulation will show the Display Form as shown in Figure 2

Display Form - Operation

Display 2 Display 2

The Display Form is divided in 5 parts

A top header section which shows general information.
A top box left showing information about 10 mini BH's
A top box right showing information about 10 mini BH's
A bottum box left showing information about 10 mini BH's
A bottum box right showing specific information about mini BH 1. This infomation is related to the first line of each of the other 3 Boxes.

Now we go to something more specific data.

The Column marked Date shows 10 dates. This shows the information the last time when the BH was either at aphelion (farthest distance) or at perihelion (closest) distance. For line 1, BH 1, this in the year 3247.
When you select Display 2 you get a different image. Now line 1, BH 1, this changes in the year 2254.
In that same column there is a section marked # of revs . This shows the total number of revolutions of the BH since the start of the simulation. The revolutions are counted based on the number of times when the mini BH reaches aphelion. For BH 1 this are 75 revolutions. When Display 2 is selected this changes to 14.
The bottom box right shows a column marked Year. Specific of interest are the 4 values in the row starting with "t 0", "t -1", "t -2" and "t max". This is the angle when the BH is at aphelion. Specific this is the longitude of the pericenter, W, which is the sum of omega and w. In this particular case for BH 1 omega = 228 degrees and w = 293 degrees (See bottum left box) which gives a total of 521 degrees or 161 degrees, which is a good match.
The reason for using w in stead of w is because the original program was used to simulate the planets around the Sun.
The column "time in sec" shows the time when the BH was aphelion. This was in cycle 193597952, 193597951 and 193597950. To calculate these numbers you have to divide the "time in sec" value by 200 which is delta time dt between each calculation cycle. This value is shown in the header behind dtime.
Using the first three values at "t 0","t -1" and "t -2" it is possible to calulate the maximum distance which is shown in the row "t max". Unfortunate the actual distances shown in column dmax don't demonstrate this.
This reflects work in progress

Display 3a

Display 3b

Both Display 3a and 3b show the same situation. The situation enlarged with a factor 2.

When Display 3a is not selected the display shows the situation in the year 2213.
When Display 3a is selected the display shows the situation in the year 16050.
When Display 3b is not selected the display shows the situation in the year 2218.
When Display 3b is selected the display shows the situation in the year 16171.

The most remarqueble chance is in the star S62, in green. Also interesting is the star S2. What is observed is the mutual influence between the 10 stars or 10 mini Black Holes.
The simulation is solely done using Newton's law.

Menu Form - Operation

Display 4

Figure Form - Operation

Display 5

VB2019 Sagittarius Software Description

Source and destination
The program supports 1 BH in the center and 10 sun's or mini BH's. Together 11 objects.
The names of the of the objects are stored in the array BHname0$(0-10)
This array stores the names in a standard order from low (S1) to high (S62)
The numbers of BH's are stored in the array source(0-10). The standard order of these numbers are from 0 to 10.
Database
The VB2019 application supports 3 databases or files. The main reason is to save the results of a specific simulation and to allow to temporarely stop a simulation and to continue at a later date.
- A Command data base. This file contains only one record. The file contains the parameters of a simulation. The two most important parameters are Subtest and test
  The two subroutines used are READCOM and WRITECOM. The filenames used end with CMn. The letter n identifies the subtest.
- A DATA file. This file contains only one record. The file contains all the parameters for all the objects, like position and velocity to continue with a simulation exactly in the same state where it was temporarily halted. The two subroutines used are READDATA and WRITEDATA
- A DataBase file. This file contains n records and is updated at regular intervals. The two subroutines used are READDB and WRITEDB
What is following is more specific information about the three files used.
- READCOM en WRITECOM use a BINARY file structure.
  At the start of each subroutine the file is opened with: FileOpen(1,filenmcm$,OpenMode.Binary.
  To read a record READCOM uses the statement FileGet(1,Comrecord
  To write a record WRITECOM uses the statement FilePUT(1,Comrecord
  At the end the file is closed with FileClose(1)
- READDATA en WRITEDATA use a Random file structure.
  At the start of each subroutine the file is opened with: FileOpen(2,filenmou$,OpenMode.Random,,,ReclenDATA).
  To read a record READDATA uses the statement FileGet(2,Datarecord, recordnumberou)
  To write a record WRITEDATA uses the statement FilePUT(2,Datarecord, recornumberou)
  At the end of each subroutine the file is not closed
- READDB en WRITEDB use a BINARY file structure.
  At the start of each subroutine the file is opened with: FileOpen(3,filenmdb$,OpenMode.Binary.
  To read a record READDB uses the statement FileGet(3,Bufrecord
  To write a record WRITEDB uses the statement FilePUT(3,Bufrecord
  At the end the file is closed with FileClose(3)
Angular orbital elements
The background is the book Explanatory supplement to the astronomical almanac paragraph 1.412 Orbital Parameters, page 27, specific Figure 1.412.1 Angular orbital elements.
In this section we use instead of mini BH the word planet. In Figure 1.412.1 the path of the planet is an ellipse and drawn counter clock wise. This defines the the orbital plane. The Blackhole is at point M. The path results in 4 states:
1. State 1 is the point N. This is point when the planet crosses the reference plane. There are two of these points (See state 3). The major characteristic of point N is that the distance towards the BH at point M (object 0) is decreasing. Point N defines the parameter omega which has a value between 0 and 360 degrees.
2. State 2 defines point P when the planet is at Perihelium or shortest distance from the BH
  Point P defines the parameter w or the argument of Pericenter, which has a value between 0 and 360 degrees.
  Point P also allows the calcultion of the inclination which has a value between -180 and 180 Point P is also used to calculate the minimum distance from the BH.
3. State 3 is the point N'. This is point when the planet crosses the reference plane. There are two of these points (See state 1). The major characteristic of point N is that the distance towards the BH at point M (object 0) is increasing. Point N' is not used to define omega.
4. State 4 defines the point when the planet reaches Aphelium or the furtest distance from the BH.
  Point P is also used to calculate the maximum distance from the BH. When both the maximum distance and the minimum distance are calculated the eccentricity can be updated.

id1 	a 	e 	i (°) 	omega  	w (°) 	Tp (yr) 	P (yr) 	Kmag 	q (AU) 	v (%c) 	dv      m0
S1 	0.5950 	0.5560 	119.14 	342.04 	122.30 	2001,800 	166.0 	14.70 	2160.7 	0.55 	0.03    12.40
S2 	0.1251 	0.8843 	133.91 	228.07 	66.25 	2018,379 	16.1 	13.95 	118.4 	2.56 	0.00	13.6
S4 	0.3570 	0.3905 	80.33 	258.84 	290.80 	1957,400 	77.0 	14.40 	1779.7 	0.57 	0.01    12.2
S6 	0.6574 	0.8400 	87.24 	85.07 	116.23 	2108,610 	192.0 	15.40 	860.3 	0.94 	0.00     9.2
S8 	0.4047 	0.8031 	74.37 	315.43 	346.70 	1983,640 	92.9 	14.50 	651.7 	1.07 	0.0     13.2
S9 	0.2724 	0.6440 	82.41 	156.60 	150.60 	1976,710 	51.3 	15.10 	793.2 	0.93 	0.02     8.2
S12 	0.2987 	0.8883 	33.56 	230.10 	317.90 	1995,590 	58.9 	15.50 	272.9 	1.69 	0.01     7.6
S13 	0.2641 	0.4250 	24.70 	74.50 	245.20 	2004,860 	49.0 	15.80 	1242.0 	0.69 	0.01    10.
S14 	0.2863 	0.9761 	100.59 	226.38 	334.59 	2000,120 	55.3 	15.70 	56.0 	3.83 	0.06    10.
S62     0.0905 	0.9760 	 72.76 	122.61 	 42.62 	2003,330 	 9.9 	16.10 	16.4 	7.03 	0.04    10.

Reflection 1

What the result of the simulations show is that the configuration is very stable.

Reflection 2

Feedback

None

Created: 8 August 2020
Updated: 15 November 2020

Back to my home page: Contents of This Document